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Introduction 

The quantum inverse scattering method was initiated some twenty years ago 
by E. K. Sklyanin and L. D. Faddeev [||, || and then mostly developed by the 
group of the Steklov mathematical institute at Leningrad (cf, for instance, 
0, §, ^, 1^). A pedagogical account of its most important aspects can be 
found in 0]. 

The quantum inverse scattering method got its name, since it arose as 
an attempt to develop a quantum version of the (classical) inverse scatter- 
ing method R M, which was successful in solving non-linear classical evolu- 



tion equations, such as the Korteweg-de-Vries equation |10|, the non-linear 



Schrodinger equation [jTI| or the sine-Gordon equation [jT2|. 

The classical inverse scattering method provides a mapping from a set of 
field variables satisfying non-linear evolution equations to a set of scattering 
data of an associated auxiliary problem. While the fields obey non-linear 
evolution equations, the scattering data evolve linearly in time. The solution 
of the initial value problem for the original non-linear evolution equations of 
the fields is achieved by first mapping the initial data to the scattering data 
at time t = 0, then using the linear time evolution of the scattering data, 
and finally applying the inverse transformation |jl^, |l4l from scattering data 



to fields at a time t > 0. 

In this paper we solve the 'inverse scattering problem' for quantum lattice 
models. The solution is remarkably simple. 

Nowadays the term 'quantum inverse scattering method' usually refers to 
a method formulated for systems of finite length. The relation to the classical 
case is the following. The elements of the monodromy matrix, which appears 
in the formulation of the classical problem for systems of finite length, have 
simple Poisson brackets [Q. In the quantum case the Poisson brackets are 
replaced by commutators of quantum operators. These commutators remain 
simple after quantization. The quantum operators can be grouped into a 
matrix, which, by analogy to the classical case, is called (quantum) mon- 
odromy matrix. The elements of the quantum monodromy matrix obey a 
set of quadratic relations. They generate the so-called Yang-Baxter alge- 
bra. The structure of this algebra is determined by numerical functions of a 
complex spectral parameter, which again can be arranged in a matrix. This 
matrix is called the i?-matrix. It satisfies the famous Yang-Baxter equation 
(see ( ^OD below). The i?-matrix and its associated Yang-Baxter algebra are 
the key concepts of the quantum inverse scattering method. These concepts 



are algebraic. 

The Yang-Baxter algebra has two primary apphcations. First of all, it 
contains, in general, a rich commutative subalgebra generated by the trace of 
the monodromy matrix. The elements of this subalgebra have a natural in- 
terpretation as a set of commuting quantities belonging to a physical system. 
One of these quantities is interpreted as the Hamiltonian. The existence of a 
large set of commuting operators can not be directly utilized to diagonalize 
the Hamiltonian. In many cases, however, the Yang-Baxter algebra can be 
employed for this task. It can be used to simultaneously diagonalize all of 
the commuting quantities by a procedure called the algebraic Bethe ansatz 
. This is the most important application of the Yang-Baxter algebra. 

Inspite of the conceptual differences between classical and quantum in- 
verse scattering method, both methods have an important point in common. 
They essentially rely on a mapping from local field variables to a set of non- 
local variables, which are the elements of the monodromy matrix. In the 
quantum case the inverse transformation, expressing the local fields in terms 
of the elements of the monodromy matrix, was not known until recently. It 
first appeared in the examples of the inhomogeneous XXX and XXZ spin-| 



Heisenberg chains in the article [|15|. 



The article ||T5[ is part of a series of articles |T6|, 0, |T^, |I8| by Izergin, 
Kitanine, Maillet, Sanches de Santos and Terras. In this series an inter- 
esting new device, the 'factorizing F-matrix' [T^ was introduced into the 
algebraic Bethe ansatz and its features were explored. This led to simplified 
derivations of a number of important results for the XXX and XXZ spin-| 
chains. Among the rederived results are the norm formulae M, 120] and the 



Slavnov formula [21] for the scalar product of a Bethe ansatz eigenstate with 
a non-eigenstate. The articles |TB|, [T^ |T^ also provide simplified derivations 
of various results for form factors [^ and correlation functions p3 
and their generalization to the case of non-zero magnetic field. 



24, 25 



From our point of view, the most interesting new result in [0, |T5|, [T^, [18 
is the solution of the quantum inverse problem for the periodic, and inho- 
mogeneous XXX and XXZ spin-| Heisenberg chains (cf [0). This result 
appeared to be the most important new tool used in the rederivation of the 
determinant formulae for form factors WE\ and in the derivation of multiple 
integral representations of correlation functions at finite magnetic field |]18| . 



In this article we shall focus on the quantum inverse problem. We shall 



obtain an explicit solution, which is valid (i) in the homogeneous caseQ, (ii) 
for models with i?- matrices of arbitrary higher dimension, and (iii), most 
generally, for fundamental graded models EH]. Upon specification to the 



cases of the inhomogeneous XXX and XXZ spin-i Heisenberg chains our 
result reduces to the formula obtained in |[T5[] . 

Our result in its most general form is given by the formula ( |5BD below. 
This formula is valid for homogeneous as well as for inhomogeneous mod- 
els. Important special cases considered in this article are the solution of the 
quantum inverse problem for the translationally invariant XYZ spin chain 
(cf equations (PTD-(P^) in section 1) and for the inhomogeneous t-J model 
(equations (|9TD-(p8D in section 6). 

Formula ( ^61) expresses the local operators as products of the entries of 
the monodromy matrix evaluated at the inhomogeneities. The structure of 
the solution of the quantum inverse problem for periodic lattice models is 
thus much different from the structure of the solution of the classical inverse 
scattering problem. In the quantum case we have an explicit multiplicative 
formula. In the classical case the solution is implicit and additive. It reduces 



to the Gelfand-Levitan-Marchenko integral equations |p!3|, P 



The article is organized as follows. In section 1 we present the solution of 
the quantum inverse problem for the homogeneous XYZ spin-| chain. Later 
this solution will appear as a special case of our general solution (^). We 
treat the case of the XYZ chain separately, since the proof greatly simplifies in 
the homogeneous case. In section 2 we remind to the reader the definitions of 
graded vector spaces and graded associative algebras. We introduce the no- 



tion of graded local projection operators which were recently defined in p6 |. 
Section 3 reviews the construction of the Yang-Baxter algebra for fundamen- 
tal graded models [^ and some important results about the graded version of 
the quantum inverse scattering method p, ^ Q. In section 4 we introduce 
canonical Fermi operators into the formalism. In section 5 we illustrate the 
abstract formalism developed in the preceding section through two examples 
which are important in physical applications. We consider the small polaron 
model |2^ and the supersymmetric t-J model |^, ^, ^, 0, |33|, pl|, p5|, |36| . 
In section 6 we present our main result (^6]) in its most general form, valid for 



inhomogeneous, fundamental graded models associated with i?-matrices of 
arbitrary dimension. We specify our formula for the examples considered in 



*The proof of |l5) does not work in the homogeneous case, since the F-matrices are 
not invertible for homogeneous lattices. 



section 5 and work out its homogeneous limit. In section 7 we give a general 
definition of the fermionic i?-operator associated with a fundamental graded 
representation of the Yang-Baxter algebra. A fermionic -R-operator was re- 

for a number of models important in physical 
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cently introduced in |^, 
applications. The fermionic i?-operator is one of the tools we shall need for 
the proof of the solution (^) of the quantum inverse problem. Section 8 is 
devoted to this proof. We construct the shift operator for inhomogeneous, 
fundamental graded models and work out its properties. Our proof of the 
quantum invers problem solely relies on the properties of the shift operator. 
In particular, we do not use factorizing F-matrices (which are so far known 
only for the XXX and XXZ spin-| chains) as in the article |jl5|. This makes 
our approach more general and powerful. Our article is concluded with a 
brief summary and a discussion of our new formulae. 



Solution of the quantum inverse problem 
for the XYZ chain 



In this section we shall start solving the quantum inverse problem for fun- 
damental models by considering an important example. The result of this 
section will later appear as a special case of our main result, equation (|86D . 
We think, however, that the structure of our main result and of its proof 
is best understood by considering an example first. We shall assume that 
the reader is familiar with the basic ideas of the quantum inverse scatter- 
ing method. Readers not familiar with those ideas are referred to section 3, 
where a brief review is provided. 



The XYZ spin-i chain is characterized by its i?-matrix |^, ^ , 



R(u) 



(a{u) d{u)\ 

h{u) c{u) 

c{u) h{u) 

\d{u) a{u)J 



(1) 



In a normalization, which assures the unitarity of the i?-matrix (cf (^)), the 



Boltzmann weights a{u), . . . , d{u) are given by 

, , sn(-u + 2r]) , , , sn(-u) 

a[u) = — — — — — , b[u^ — 



sn{u) + sn{2ri) ' sn('u) + sn{2ri) 

(2) 
sn{2ri) k sn{u)sn{2ri)sn{u + 2ri) 

du) = — 1-^ — r , dlu) = — ^ — r 

sn[uj + sn[2ri) sn{U) + sn(2?7) 

The i?- matrix R{u) is considered as acting on the tensor product C^ ® C^. 
R{u,v) := R{u — v) is a solution of the Yang-Baxter equation (^Dp (see 
below). Hence, an exactly solvable spin chain can be associated with R{u). 
Let us briefly recall the steps necessary for its construction. 
Deflne e^ G End(C2), a, /? = 1, 2 by 






The set {e^ G End(C2)|a,/5 = 1,2} is a basis of End(C2). The deflnition 

e,f = Jf ^-^) ® e^ ® Jf ^-) , (4) 

for j = 1, . . . ,L, and I2 being the 2x2 unit matrix, provides a basis of 
End(C^)'^'^, which is the space of states of an L-site spin-| quantum spin 
chain. The matrices e^^ satisfy 

[e,f,e4] = ioTj^k , e,fe4 = 5?e,^ (5) 

and have the meaning of local projection operators. Using the i?-matrix (|l]) 
and the local projection operators we can deflne the L-matrix at site j as 

(6) 

The L-matrix Lj{u) is a 2 x 2 matrix in an auxiliary space. Its entries 
are operators acting on the space of states of an L-site spin-i chain. The 



monodromy matrix of the corresponding homogeneous spin chain is the L- 
fold ordered product 

T{u) = Ll{u)...Li{u) . (7) 

By construction the monodromy matrix gives a representation of the Yang- 
Baxter algebra with i?-matrix R{u) (see section 3 below). We have 

PR{u - v) {T{u) ® T{v)) = {T{v) ® T{u))PR{u - v) , (8) 

where P = Y2^ g^^ e^ ® e^ is the permutation matrix on C^ (S> C^. The Yang- 
Baxter algebra (^ is the basis for the solution of the XYZ chain by algebraic 
Bethe ansatz 0. 

The monodromy matrix (|^) is a 2 x 2 matrix in auxiliary space. Its entries 
are non-local operators acting on the space of states of the XYZ chain. We 
may write the monodromy matrix as 

[a(u) Biu) , 
\C{u) D{u) ' 

The quantum inverse problem is to express the local operators e„^ in terms 
of the elements A{u), . . . , D{u) of the monodromy matrix. 

For the case at hand this problem is rather easily solved. Note that 
a(0) = c(0) = 1 and 6(0) = d{0) = 0. It follows from (D that 

^.(0)= ('^1 'i] = Y.ei®er, = Po, . (10) 

Here Poj is the permutation operator that interchanges the auxiliary space 
with the j-th quantum space. Similarly, Pjk = J2a /3=i ^ia^fc/3 interchanges 
the j-th and k-th quantum spaces. Using equation (|1^) in equation (|^) we 
obtain 

T{0) = PolPol-i...Poi=PoiPilPil-i...Pi2 = PoiU . (11) 

Here U = PilPil-i ■ ■ ■ P12 = P12P23 ■ ■ ■ Pl-il is the cyclic shift operator in 
quantum space. In the second equation in (11) we have used the commutation 



relation PjkPji = PkiPjk for permutation operators. Let us write equation 
([TT|) in matrix form, 



(12) 



A{0) 5(0) \ ^ lei\U ei\U 
C{0) D{0)) ~ [e,lU e,lU^ 
Comparing the matrix elements we find the relations 

U = A{0) + D{0) , (13) 

a- = e,l = B{0)U-' , (14) 

at = e,j = CiO)U-' , (15) 

at = e,\-eil = {A{0)-D{0))U-' . (16) 

These relations constitute a solution of the quantum inverse problem for 
the local operators acting on the first lattice site. We may now simply use 
the shift operator to shift the site indices. Since U"'~^ei^U^~"' = e„^ and 
U^ = id, we obtain 

a- = f/"-ifi(0) f/^-" , (17) 

a+ = ?7"-^C(0) t/^~" , (18) 

a^ = f/"-i (A(0) - D(0)) U^-" . (19) 

Taking into account equation (0) we see that the right hand side of the 
equations (|T^)-([TP|) are entirely expressed in terms of the entries of the mon- 
odromy matrix. An alternative way of writing (p!^)-([T9|) is 

ej^ = t/"-i r„^(0) f/^-" = (A(0) + D{0)y-' Tf (0) (A(0) + D(O))''-" . 

(20) 

Equation ( pOD allows us to calculate expectation values of local operators by 
means of the Yang-Baxter algebra. 

The remainder of this article will be devoted to the generalization of equa- 
tion (|20|) to (i) arbitrary dimension of the i?-matrix, (ii) the inhomogeneous 
case, and (iii) to fundamental graded models. It is important noting that 
our above solution of the quantum inverse problem does not depend on the 
specific features of the XYZ chain. Our calculation solely relied on the fact 
that the L-matrix evaluated at u = turns into a permutation operator (cf 
equation (p!0D ). 



2 Graded vector spaces 

In this section we shall recall the basic concepts of graded vector spaces and 
graded associative algebras. In the context of the quantum inverse scatter- 



ing method these concepts were first used by Kulish and Sklyanin g, ^ . 
We shall further recall the notions of 'graded local projection operators' and 
graded permutation operators. Graded local projection operators were intro- 
duced in the article [^. They enable the definition of fundamental graded 



representations of the Yang-Baxter algebra, which will be given in the fol- 
lowing section. 

Graded vector spaces are vector spaces equipped with a notion of odd and 
even, that allows us to treat fermions within the formalism of the quantum 
inverse scattering method. Let us start with a finite dimensional local space 
of states V, on which we impose an additional structure, the parity, from 
the outset. Let V = Vq Q) Vi, dimVo = m, dimVi = n. We shall call 
fo £ K) even and f i G Vi odd. The subspaces Vq and Vi are called the 
homogeneous components of V. The parity p is a function V^ — > Z2 defined 
on the homogeneous components of V, 

p{vi) =i , i = 0,l , VieVi . (21) 

The vector space V endowed with this structure is called a graded vector 
space or super space. Let us fix a basis {ei, . . . , Cm+n} of definite parity and 
let us define p{a) := p{ea)- 

The use of graded vector spaces within the quantum inverse scattering 
method requires the construction of an algebra of commuting and anticom- 
muting operators. For this purpose we have to extend the concept of par- 
ity to operators in End(V^) and to tensor products of these operators. Let 
e^ G End(V"), e^e^ = S^e^. The set {e^ G End{V)\a,P = 1, . . .m + n} is a 
basis of End(l^). Hence, the definition 

p{e'i)=p{a)+p{(3) (22) 

induces a grading on End(y) regarded as a vector space. 

It is easy to see that an element A = A'aC^ G End(K) is homogeneous 
with parity p{A), if and only if 
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The latter equation implies for two homogeneous elements A,B& End(l^) 
that their product AB is homogeneous with parity 

p{AB) = p{A) + p{B) . (24) 

In other words, multiplication of matrices in End(y) preserves homogeneity, 
and, therefore, End(\^) endowed with the grading (p2|) is a graded associative 
algebra 0. 

Let us consider the L-fold tensorial power (End(V"))®^ of End(V^). The 
definition (|22|) has a natural extension to (End(V"))'^'^, namely, 

p{e1\ (» ■ ■ ■ (^ e^^^J = p{ai) + p{Pi) + ■ ■ ■ + p{aL) + pWl) • (25) 

From this formula it can be seen in a similar way as before, that homogeneous 
elements A = A^';;;|^e^; ® ■ ■ ■ ® e^^ of (End(\/))®-^ with parity p{A) are 
characterized by the equation 

(_l)Et=i(pK)+p(/3.))A;^i---^i = (-l)P(^)A^i;;;|^^ , (26) 

which generalizes (|2^). Again the product AB is homogeneous with parity 



p{AB) = p{A) +p{B), if A and B are homogeneous. Thus the definition (pSf ) 
induces the structure of a graded associative algebra on (End(l^))®^. 
Let us define the super-bracket 

[X, y]± = XY- (-i)pwp(^)rx (27) 

for X, Y taken from the homogeneous components of End(l^), and let us 
extend it linearly to End(V^) in both of its arguments. Then, End(y ) endowed 
with the super-bracket becomes the Lie-super algebra g\{m\n). Note that the 
above definition of a super-bracket makes sense in any graded algebra and is 
particularly valid in (End(V^))®^. 

The following definition of 'graded local projection operators' [^ will 



be crucial for our definition of fundamental graded representations of the 
Yang-Baxter algebra in the next section. Define the matrices 

e^./3 = (_i)(pW+p(/5))EL.+ip(7.) jg;^) ® e^ ® e:^ ® ■ ■ ■ ® e^^ , (28) 

where Im+n is the {m + n) x [171 + n) unit matrix, and summation over double 
tensor indices (i.e. over 7j+i, . . . ,7l) is understood. We shall keep this sum 
convention throughout the remainder of this article. The index j on the left 
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hand side of ( pS]) will later refer to the j-th site of a physical lattice model 
and is called the site index. A simple consequence of the definition (pSD for 
j y^ k are the commutation relations 



e.fe4 = (-l)(^(")+^(^«(^(^)+^(^»e4e,f . (29) 

It further follows from equation (123) that Cj^ is homogeneous with parity 

p{e/j = pia) + piP) . (30) 

Hence, in agreement with intuition, equation (123) says that odd matrices with 



different site indices mutually anticommute, whereas even matrices commute 
with each other as well as with the odd matrices. For products of matrices 
Cj^ which are acting on the same site ( [28| ) implies the projection property 



P^S _ A/3..5 



e,„e 



Ja^J-r "7 



<e,^ . (31) 



Equations (^) and (^Tj) justify our terminology. The Cj^^ are graded ana- 
logues of local projection operators. We call them graded local projection 
operators or projection operators, for short. Using the super-bracket (P^, 
equations (pOj) and (|31|) can be combined into 



[e.f , e4]± = 6,, (6^^e,i - (-l)M'^)^p(mvi.)^PiS)) ^s^^A . (32) 



The right hand side of the latter equation with j = k gives the structure 
constants of the Lie super algebra gl(m|n) with respect to the basis {cj^}. 

Since any m + n-dimensional vector space over the complex numbers is 
isomorphic to C^^", we may simply set V = C"^'^"'. We may further assume 
that our homogeneous basis {e^ G C"^~^"'\a = 1, . . . ,m + n} is canonical, i.e. 
we may represent the vector Cq by a column vector having the only non-zero 
entry -fl in row a. Our basic matrices e(^ are then {m + n) x (m + n)-matrices 
with a single non-zero entry +1 in row a and column j3. 

Remark. The meaning of (p8[) becomes more evident by considering a sim- 
ple example. Let m = n = 1 and p(l) = 0, p(2) = 1. Then, using (|5^), we 
obtain 

[cji, efci]± = {cjl, Cfci) = , (33) 

[ej2' ^kl]± = {ejl, Ckl} = , (34) 
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[Gji^ Gk2\± — |Cj25 6fci| — Ojk{eji + ejg) — Ojk (35) 



for j, fc = 1,...,L. The curly brackets in (p3D, (0) denote the anticom- 
mutator. The matrices Cj^ and ek\ satisfy the canonical anticommutation 



2 



C 



-J 



relations for spinless Fermi operators. We can therefore identify Cj^ 
and efc2 — ^ c^- Introducing Pauh matrices a^ = ef, o~ = e\ and a^ = ^\~^2 
we obtain, by carrying out the summation, the following explicit matrix rep- 
resentation from our basic definition (|28|), 



c 



'3 



jf •'■-'' ® a+ ® (a^)®(^-^') , (36) 



c+ = /f'='i)®a-®(a^)®(^-'=) . (37) 

This is the well-known Jordan- Wigner transformation |^1[] expressing Fermi 
operators for spinless fermions in terms of Pauli matrices. We may thus inter- 
pret equation (^) as a generalization of the Jordan- Wigner transformation. 
In general, equation ( pSD provides matrix representations not of Fermi oper- 
ators but, more generally, of fermionic projection operators. Representations 
of Fermi operators can be obtained be taking appropriate linear combinations 
of matrices e^^. This point will be elaborated in section 4 below. Note that an 
alternative generalization of the Jordan- Wigner transformation, which relies 



on representations of the Clifford algebra, was recently introduced in [Q . 

The permutation operator plays an important role in the construction of 
local integrable lattice models. It enters the expression for the shift operator 
on homogeneous lattices and the expression for the Hamiltonian of rational 
gl(m|n) invariant models (see section 5 below). In the graded case the def- 
inition of the permutation operator requires the following modifications of 
signs, 

P,. = (-ir^^^e.fe,^ . (38) 

As indicated by its name, this operator induces the action of the symmetric 
group ©^ on the site indices of the matrices Cj^ . The properties of Pjk (for 
j 7^ k) are the same as in the non-graded case. They are easily derived from 
(|29|) and ( plj) and can be found, for instance, in pGf. Let L = 2. Then 



(39) 

From the right hand side of this equation we can read off the matrix elements 
of Pi2 with respect to the canonical basis of End(l^ V). 
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Figure 1: The Yang-Baxter equation is most easily memorized in graphical 
form 

3 Fundamental graded models 



In this section we shall recall the notion of fundamental graded representa- 
tions of the Yang-Baxter algebra, which was recently introduced in p6[. For 



a given grading we shall associate a fundamental model with every solution 
of the Yang-Baxter equation that satisfies a certain compatibility condition 
(see @ below). 

For our present purpose it is most suitable to interpret the Yang-Baxter 
equation as a set of functional equations for the matrix elements of an {m + 
n)^ X [m-\- ra)^-matrix R{u, v). We may represent the Yang-Baxter equation 
in graphical form as shown in Figure 1, where each vertex corresponds to a 
factor in the equation 



^S' (^' ^)^a"7' (^' ^)^/3"7" I 



V, w) 



Rf}7^' ( 



^)K''y" 



(n,u;)i?°„^„(n,w) 



(40) 



Note that there is a direction assigned to every line in figure 1, which is 
indicated by the tips of the arrows. Therefore every vertex has an orientation, 
and vertices and i?-matrices can be identified according to figure 2, where 
indices have been supplied to a vertex. Summation is over all inner lines in 
figure 1. 

The construction of a graded Yang-Baxter algebra and its fundamental 
representation requires only minimal modifications compared to the non- 
graded case 0. Let us assume we are given a solution of (|40|) , which is 
compatible with the grading in the sense that ||^ 



<(«,^) 



_l\p{a)+p{l3)+p{t)+p{5) j^af3 ^^ ^ 



(41) 
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<(«,^) 



u 



Figure 2: Identification of tlie i?-niatrix witli a vertex 
Define a graded L-matrix at site j as 



(42) 



Equation (^If ) implies that the matrix elements of !ij{u,v) are of definite 
parity, 

p{L,^^{u,v))=p{a)+p{(3) . (43) 

Thus their commutation rules are given by 

Si^l{u,v)L,]{w,z) = (-l)(P(°)+^(^))(^(^)+^(^))£fc](t/;,z)£,;(w,t;) . (44) 



It further follows from the Yang-Baxter equation (^) and from equation (^if) 
that 

R{u,v)(^fij{u,w) ®s ^j{v,w)) = (^Cjj{v,w) ^s iij{u,w))R{u,v) . (45) 
As in the non-graded case the matrix R{u, v) is defined by 

<K^) = <K^) • (46) 

The super tensor product f^ in equation (|45|) is to be understood as a su- 
per tensor product of matrices with non-commuting entries, {A (S)^ B)'^J = 
(^_l)(p(")+p(/3))p(7)^a2jj^ The super tensor product has the following impor- 
tant feature. Given matrices A, B, C, D with operator valued entries, which 
mutually commute according to the same rule as Lj and L^m equation 
we obtain for the product of two super tensor products 



{A ®s B){C 0s D) = AC 0s BD 



(47) 
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Equation (^Sf ) may be interpreted as defining a graded Yang-Baxter algebra 
with i?-niatrix R. We call iLj its fundamental graded representation. 

Starting from ( ^5]) we can construct integrable lattice models as in the 
non-graded case 0. Let us briefly recall the construction with emphasis 
on the modifications that appear due to the grading. Define a monodromy 
matrix 7{u, v) as an L-fold ordered product of fundamental L-matrices, 

7{u,v) = liL{u,v)...Li{u,v) . (48) 

Due to equation (^) the matrix elements of 7{u, v) are homogeneous with 
parity p(T^(m, v)) = p{a)+p{(3). Repeated application of p5| ) and (^) shows 
that this monodromy matrix is a representation of the graded Yang-Baxter 
algebra, 

R{u, v) (T(u, w) ®s 7{v, w)) = {7{v, w) (g)^ 7{u, w))R{u, v) . (49) 

In the non-graded case (n = 0) the super tensor product in (|4^) agrees with 
the usual tensor product. Let us now define the super trace as 

str(v4) = (-1)P("M;^ . (50) 

It follows from (|1|) and (||) that 

[sti{7{u,w)),stii7{v,w))]=0 , (51) 

which is in complete analogy with the non-graded case. 

Let us assume that R{u, v) is a regular solution of the Yang-Baxter equa- 
tion, i?"f (t;, v) = 626^. Then (H) implies that 

£,^(t;,t;) = (-l)^(")^('3)e,; , (52) 

and we can easily see [EH] that the super trace of the monodromy matrix 



evaluated at (v, v) generates a shift by one site, 

StT{7iv,v)) = Pi2P23.--PL-lL=:U . (53) 

It follows that t{u) := — iln(str(T(M,f ))) generates a sequence of local oper- 
ators [43] which, as a consequence of (|5TD , mutually commute. 



T{u) = tl+{u-v)H + 0{{u-vf) . (54) 
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n in this expansion is the momentum operator. On a lattice, where the mini- 
mal possible shift is by one site, and thus U rather than 11 is the fundamental 
geometrical operator, some care is required in the definition of 11. As was 



shown in [Q a proper definition may be obtained by setting 11 := — iln([/) 
mod 27r and expressing the function f[x) = x mod 27r by its Fourier sum. 
Then 11 becomes a polynomial in U. 



m=l \ / 

where = 2tt/L. The first order term H in the expansion ( p^ may be 



interpreted as Hamiltonian. Using (^) it is obtained as 



L 



H = J2Hjj+i , (56) 

where -f^LL+i = Hn and 

iJ,,+i = -i (-1)^(7)(P(«)+P(7)) 0^ R-J{u, v) _ e,^e,+4 . (57) 



u=v 



We would like to draw the readers' attention to the following points, 
(i) The i?-matrix R in equation (|^) does not undergo a modification due 
to the grading, (ii) The only necessary compatibility condition which has 
to be satisfied in order to introduce a fundamental graded representation 
of the Yang-Baxter algebra associated with a solution of the Yang-Baxter 
equation is equation (^T]), which was introduced in 0. Equation ( PD is a 
rather weak constraint. The most important solutions of the Yang-Baxter 
equation, which appear in physical applications are compatible with a non- 
trivial grading (see the examples in section 5). Moreover, a given i?-matrix 
may be compatible with different gradings, leading to different fundamental 
graded representations of the Yang-Baxter algebra ||26|| . 

Before turning to our next subject let us introduce the inhomogeneous 
generalization 

T(n;ei, . . . ,a) = ^l{u,^l) ■ ..^i{u,^i) (58) 



of the monodromy matrix (^§1). This monodromy matrix satisfies (^). For 
^i = ■ ■ ■ = C,L = V it turns into 7{u, v) defined in (|i8|). We shall formulate our 
main result below for the inhomogeneous model generated by T(u; ^i, . . . , C,l)- 
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4 Fermi operators 



In the article [^ it was explained how the various graded objects introduced 
in the previous section can be expressed in terms of Fermi operators. The key 
observation is that, as far as the matrices e,-^ are concerned, all calculations 



of the previous section rely on the commutation relations (|29D and on the 



projection property (31). Fermionic projection operators satisfy the same 



equations. We may thus say that the matrices e^^ are matrix representations 
of fermionic projection operators. As we have seen in the previous section, 
the matrices Cj^ are suitable for formulating a graded version of the quantum 
inverse scattering method. For the physical interpretation of the models 
constructed from a given solution of the Yang-Baxter equation, however, it 
is convenient to introduce Fermi operators into the formalism. 

A general construction of fermionic projection operators for fermions of 
arbitrary su(N) spin was presented in [^. Rather than repeating those 
results let us illustrate them by example. 

Let us consider spinless fermions on a ring of L lattice sites, 

{cj,Cfc} = {c],4} = , {cj,cl}=6jk , j,k = l,...,L . (59) 
It is easy to verify that the entries (Xj)'^ of the matrix 




(60) 

are fermionic projection operators. Define Xj^ = {Xj)g. Then 

X/X/ = S^X/ . (61) 

The operators Xj^ carry parity, induced by the anti commutation rule (^) for 
the Fermi operators. For j ^ k Xj^ and X^ anticommute, if both are build 
up of an odd number of Fermi operators, and otherwise commute. This fact 
can be expressed as follows. Let p{l) = 0, p(2) = 1 and p{Xj^) = p{a)+p{/3). 
Then Xj^ is odd (contains an odd number of Fermi operators), iip(Xj^) = 1, 
and even, if p(Xj^) = 0. The commutation rules for the projectors Xj are 
thus 

X/^Xkt, = (-i)(p(")+p(/3))(p(7)+PW)Xfc^X,-;^ . (62) 
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respectively. Since 
j) and (0), we may 



Now (0) and (^) are of the same form as (|3T|) and 
the calculations in the previous section relied solely on 
simply replace Cj^ — >• X^-^ in equations ( ^2] ) and (|57D. 

Fermionic representations compatible with arbitrary grading can be con- 
structed by considering several species of fermions and graded products of 
projection operators. We shall explain this for the case of two species. This 
is the case most interesting for applications, since we may interpret the two 
species as up- and down-spin electrons. We have to attach a spin index to the 

Co 



Fermi operators, Cj — >• Cjo-, a =],i, {cja^cl^} = ^jk^crr- Accordingly, there 
are two species of projection operators, Xj^ -^ -^j ■ 

Let us define projection operators for electrons by the tensor products 



J a-y 



_^N(p(a)+p(/3))p(7)j^i/3^T'5 



X 



x\ 



07 



p& 



Then 



Y 0^ Y ^'^' - Ji^ Ji^ Y P'^' 



(63) 



(64) 



Xj^ inherits the parity from X and X- . The number of Fermi operators 



contained in Xj^ is the sum of the number of Fermi operators in X|^ and 









Xj^. Hence p{X/J^) = p{XjJ +p(XJ^) = p{a) + ■ ■ ■+p{S), and the analog 



of (|62D holds for Xj^ , too. Again we present all projection operators in form 



of a matrix (X 



07 
07 X -^^ 

3 a'y 



jjps 



X, 



x\ ®. X] 



(1 - nj{)c\^ 



(1 - nj^)nj^ -Cjic]^ 



\ 















n^ 



(65) 



i%T / 



Here we used the standard ordering of matrix elements of tensor products, 
corresponding to a renumbering (11) -^ 1, (12) -^ 2, (21) -^ 3, (22) -^ 4. 
Within this convention Xj^ is replaced by Xj^^ a, /5 = 1, . . . , 4, which then 



satisfies (|6l| ) and (|62D with grading p(l) = p(4) = 0, p{2) = p(3) = 1. 



19 



Note that Fermi operators can be obtained as linear combinations of 
projection operators. We have, for instance, Cj-^ = Xj\ + Xj'^. 

So far we have considered the case of spinless fermions with two-dimen- 
sional local space of states and grading m = n = 1, and the case of electrons 
with four-dimensional space of states and grading m = n = 2. There are four 
different possibilities to realize (p9D and (|3TD in case of a three-dimensional 
local space of states, m + n = 3. They can be obtained by deleting the a's 
row and column of the matrix Xj in equation (|65|), a = 1,2,3,4. (|6TD and 



(|6^) remain valid, since the operators Xj'^ are projectors. 

An alternative way [^, ^, ^ of introducing Fermi operators into the 
quantum inverse scattering method is by applying the Jordan- Wigner trans- 
formation PT| to the non-graded L-matrix and then pulling out the non- 
local factors. This approach was of primary importance, for instance, for a 
fermionic formulation of the Yang-Baxter algebra of the Hubbard model ^E 
and led to the discovery of a S0(4)-invariant form of the monodromy matrix 
of the Hubbard model [^, 0, ^]. In general, however, we prefer the method 
presented above, since the approach of ^5|, ^, ^ so far did not lead to gen- 
eral formulae such as 



or (^) and may have unpleasant side effects, such 
as a twist of boundary conditions or the appearance of numerous factors of 
'i' in the equations. 



5 Examples 



Before turning to our main result let us present several examples in order to 
provide an idea to the reader which applications we have in mind. 

We shall start with the XXZ spin chain. The i?-matrix of the XXZ spin 
chain can, for instance, be written as 



R{u,v) 



/l o\ 

b{u,v) c{u,v) 
c{u,v) b{u,v) 



\0 



where 



b{u,v) 







sh.{u — v) 
sh{u — V + in) 



, c{u,v) 



(66) 



V 



sh(ifi;) 
sh('U — V + in) 



(67) 



20 



The i?-matrix (^) is compatible with the grading p{l) = 0, p{2) = 1 (cf 
equation (^I]) ). The corresponding L-matrix then follows from (^21), 



ejl + b{u,v)ejl c{u,v)ejl 
Lj[u,v) =1 I . (6? 

c{u, v)ej\ b{u, v)ej\ — e,^ 



^J2 

Using (|56[) , ([57|) we obtain the Hamiltonian 



-1 ^ 
-^ = • / \ Z^|^i2^i+ii + ^i+i2^ii ~ cos(k) [cj^Cj+i^ + ej+i^ejg) | • 

(69) 
We may now replace the matrices Cj^ by the fermionic projectors Xj^ ^ equa- 



tion (|60|). Then 



c(m, f ) Cj b{u, f )(1 — rij) — rij 



and 



-1 ^ 
i^ = . , , y^{c|c,+i + cf,,Cj + 2 cosffi;) n,ra,+i| + 2ctg(K)iV , (71) 

where A^ = X]i=i^j is the particle number operator. The Hamiltonian (fflD 
defines the 'small polaron model' [^. Note that the algebraic Bethe ansatz 
for the small polaron model gets slightly modified compared to the 'non- 
graded' spin chain case (p(l) = p(2) = 0), since due to the grading p(l) = 0, 
p(2) = 1 there appear certain minus signs in the Yang-Baxter algebra (|15|). 
Our next example is the well known family p^ of graded rational R- 
matrices 

R';^{u, v) = a{u, t;)(-l)p(-)p(/3) 6^6^s + ^(«, ^) ^s^^, , (72) 

where 

u — V i 
a(u,v) = , d(u,v) = . (73) 

R{u, v) solves the Yang-Baxter equation (|40|) for arbitrary matrix dimension 
_/V^ X A^^ and arbitrary grading p : {1, . . . , N} -^ Z2. 
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Remark. Note the following subtlety. The grading introduced in (|T2D is 
independent of the grading that enters the definition (|28| ) of the matrices Cj^. 
Let q : {1, . . . , N} -^ Z2 arbitrary. Then, because of the Kronecker deltas in 

m 



<(«,^) 



'_]_^g(a)+g(/3)+g(7)+g(5) _^"/3('^^ ^ 



(74) 



i.e. the compatibility condition (|4l|) is satisfied for arbitrary p and q. For 
example, let N = 2, p{l) = p{2) = 0. Then R{u,v) is the i?-matrix of 
the XXX spin-i Heisenberg chain, which is compatible with the grading 
g(l) = 0, g(2) = 1 leading to a special case of the small polaron Hamiltonian 
introduced above. 



Let us now elaborate on the case p = q. Since a{v, v) = and d{v, v) = 1, 
the i?-matrix defined in equation (|7^ ) is regular. Furthermore, 



duR'^^{u,v] 



^"'^f 



.]_-\P(a)p{/3) ^a^P 



(75) 



Thus, using p = q in equations ( |5B| ) and (0) we find the Hamiltonian 



H 



/ .(-Pji+i 



-1^ 



(76) 



where -P,j+i is the graded permutation operator defined in (^8]). Clearly H 
commutes with the generators 



E? 



E 



'J a 



(77) 



of the graded Lie algebra gl(m|n). 

The family of Hamiltonians (|7^) based on graded permutations includes a 
number of models that are interesting for applications in physics. In the 'non- 
graded' case (p(a) = 0, a = 1, . . . , A^) we have the XXX spin-| chain and its 
su(A^) generalizations. The case m = 1, n = 2 leads us to the supersymmetric 



t-J model [^. In order to see this we shall employ the fermionization scheme 
of the previous section. We start with the set of projectors obtained from 
the matrix Xj in ( |65|) by deleting its fourth row and column. We obtain the 
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reduced matrix 



Al - ^ii)(l - ^it) (1 - nji)cjr Cjiil - nj^)\ 



X, 



(1 - ^ii)ST 



:i-n 



jD'^n 



"'^iiS'T 



(7J 



sViT nji{l-nj^)J 



The entries (^j)^ of this matrix form a complete set of projection operators 
on the local space of states spanned by the basis states |0), c^||0), c^j|0). 
Double occupancy is excluded on this space of states. Let Xj^ = (Xj)'^. The 
operator 



J a 



%T^ii 



(79) 



projects the local Hilbert space of electrons onto the space with no double 
occupancy. The global projection operator for a chain of L sites is given by 
the product 



A 






%T^iiJ 



(80) 



The permutation operator Pj^ is given by equation ( pH]) with X^^ replacing 
Cn^ . The summation in (|3^) is now over three values, a, /? = 1, 2, 3, and the 
grading is p(l) = 0, p(2) = p(3) = 1. An elegant way of taking into account 
the simplifications arising from the restriction to the Hilbert space with with 
no double occupancy is to consider Pjk^ instead of Pj^. Since rij^riji^ = 0, 
we obtain 

(Pjk - 1)A = A{c]^Cka + cLcja)A - 2{S^S^ - \njnk)A - (n, + nk)A . 



Here we have introduced the electron density rij 
densities 



S^ 



2'^af3'^ja'^jf3 



n,jf + n.ji and the spin 



^2) 



The a"", a = x,y, z, are the Pauli matrices, and we identify 1 with | an 2 with 
J, in the summation over a and jS. Inserting (|8T|) into the expression (|7^) for 
the Hamiltonian we obtain the familiar Hamiltonian of the supersymmetric 



t-J model 29, 30, 32, 33, M 
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Let us also write down the corresponding L-matrix, which follows from 

^32 



equation (^^ 



(xA 



Lj{u,v) 



a{u,v) + d{u,v) 



X 



'Xj2 



^33 * 
-Xjl 



(83) 



V^ii "^32 ~XnJ 



This form of the L-matrix suggests a similar form for the monodromy matrix 
of the corresponding inhomogeneous model, 



7{u-i,,...,iL) 



( A{u) Biiu) B2{u)\ 

Ciiu) D\{u) Dliu) 
\C,{U) Dliu) Dliu)) 



U) 



6 Solution of the quantum inverse problem 

We are now ready to formulate our main result, which is a formula that 
expresses the local projection matrices Cj^ for fundamental graded models 
in terms of the elements of the monodromy matrix. We shall assume we 
are given a solution of the Yang-Baxter equation (|40|) which is regular and 
unitary. Unitarity means that -R(m, v) satisfies the equation 



ja/?/ 



-,(57 



R;^i^,v)R"JA^,^) = sp'a' 



(85) 



Let p be a grading that is compatible with the i?-matrix in the sense of 
equation (|4l|) , and let 'J(m) = T(u; ^i, . . . , C,l) be the corresponding inhomo- 
geneous monodromy matrix (pHI). Then we have the following formula. 



n—l 



--na 



.i)PWp(/3)-Qgtr(T(e,))-o'f(en)- n ^tr(^(^.)) 



(86) 



i=i 



j=n+l 



Equation (^) is our main result. It constitutes a solution of the quantum 
inverse problem for fundamental graded models. We shall prove it in the 
remaining sections of this article. For m = 2, ra = (p(l) = p(2) = 0) 
equation (|86D reduces to a result recently obtained by Kitanine et al. ||15| . 
Note that because of ( |5lD no ordering is required for the products on the 
right hand side of (pGl). 
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Before proceeding with the proof of (|86D let us illustrate the equation 
through the examples of the previous section. Note that the functions b{u, v), 
c{u, v) in ( p7| ) and a{u, v), d{u, v) in ( [73| ) have been chosen in such a way that 
the corresponding i?- matrices (|66|), (|7^) satisfy the unitarity condition (^S]). 

For the small polaron model the monodromy matrix is of the form 

'a(u) B(u) , 
7iu) =\ ^ ^ \ ■ (87) 

\C{u) D{u)' 

Using the fermionization ( pDj ) we obtain from 

n-l L 

c+ = n(^(^^) - ^(^^■)) ■ ^(^") ■ n (^(^^■) - ^(^^•)) ' (88) 

j=l j=n+l 

n-l L 

cn = n(^(^^) - ^(^^■)) ■ ^(^") • n i^^^^) - ^(^^■)) ' (89) 

j=l j=ra+l 

n-l L 

nn = - n(^(^.) - D{Q) ■ D{i^) . n (^(e,) - D{Q) . (90) 

j=l j=n+l 

A similar set of equations holds for the local operators of the super- 
symmetric t-J model. The monodromy matrix was presented in equation 
(|8^). Since the grading is p(l) = 0, p{2) = p(3) = 1, we have str (T(-u)) = 
A{u) -tr(D(u)), and 

n-l L 

{l-nni)c+^=l[stT{7{Q).B,{U. H stT{7{Q) , (91) 

j=l j=n+l 

n-l L 

c+(l-n„t)=nstr(T(e,))-52(en)- n ^^^(^(^^O) ' (92) 

j=l j=n+l 

n-l L 

{l-nni)cn^=l[stT{7{Q)-C^i^n)- U StT{7{Q) , (93) 

j=l j=n+l 

n-l L 

Cni{l-nn^)=l[stT{7{Q)-C2{U- H StT{7{Q) , (94) 

j=l j=n+l 
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n-l L 

St = -\[M7{^i)) ■ Dli^n) ■ W M^Q) , (95) 

j=l j=n+l 

n-l L 

S- = -l[ str(0-(O)) ■ DliU . H str(T(0)) , (96) 

j=l j=n+l 



^ n— 1 L 

^" = -9 n str(T(0)) ■ (/^K^n) - /^^'(en)) ■ n str(T(0)) 



j=l j=n+l 

(97) 

n-l L 

(l-n„J(l-n„t)=nstr(7(e.))-^(U- n ^t^(^(^^)) • (98) 

j=l j=n+l 

On the restricted Hilbert space of the supersymmetric t-J model, where dou- 
ble occupancy of lattice sites is excluded, the operators on the left hand side 
of equations (|9TD-(p4D are appropriate creation and annihilation operators. 
The operator (1 — ?T-ni)c^^, for instance, creates an up-spin electron on site n, 
provided this site is not occupied by a down-spin electron. The local spin op- 
erators S^ (a = X, y, z) were introduced in equation (^2]). S^ and S~ in (|95|) 
and (|6|) are defined as 8^ = 8^ + i8l = c^^Cni and 5" = 5"^ - i8l = c+j^c„|. 
These operators induce a local spin flip. The operator (1 — ri„j)(l — n„|) on 
the left hand side of ( ]98| ) acts as 1 — (n„| -|- n„j) on the restricted Hilbert 
space of the supersymmetric t-J model and thus essentially gives the local 
particle number operator. 

Note that equations (|9lD-(|98D after appropriate replacement of the mon- 



odromy matrix also apply to the 't-0 model' (the infinite coupling limit of 
the Hubbard model below half-filling) which was solved by algebraic Bethe 



ansatz in ||26| 



We would like to stress, that our main result, equation (pq), also holds 
for homogeneous, translationally invariant models, for which C,j = v = 0, 
for j = 1, . . . ,L. In this case ( pB] ) takes a particularly simple form, since 
str(T(0)) = U, where U is the homogeneous shift operator (|53|) . Using this 
fact equation (|8^) turns into 



a 



J^ = (_i)p(")p{/3) t/«-iT^(o) ?7^-« . (99) 



Performing a similarity transformation with U " we obtain the amazingly 
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simple result 

ei^ = (-l)p(-M/3) Tf (0) f/-i , (100) 

which in section 1 was obtained for the special case of the XYZ chain. To give 
another example, equations (|88D -(pO|), for instance, are in the homogeneous 
case equivalent to 

ct = B{0)U-^ , ci = C{0)U-^ , ni = -D{0)U-^ . (101) 

7 The fermionic i^-operator 

The role of the matrix R{u, v) in the graded Yang-Baxter algebra ( ^5]) is to 
switch the order of the two auxiliary spaces. The definition of an operator 
that similarly switches the order of quantum spaces in a product of two L- 
matrices requires appropriate use of the grading. Recently, such an operator 
was introduced for several important models by Umeno, Shiroishi and Wadati 



37, 3S] and was called fermionic /^-operator. Here we give a general definition 



of the fermionic i?-operator associated with a solution R{u, v) of the Yang- 
Baxter equation (^0]) . For a given grading and a solution R{u,v) of the 
Yang-Baxter equation ( |ID| ) that is compatible with this grading (cf (^)) we 
define 

3^Jfc(«, v) = (-l)p(7)+p(")b(/3)+p(7)) R-^^u^ ^)e^.7e,J . (102) 

The fermionic i?-operator will be an important tool in the proof of our main 
result (p6[). Let us summarize its properties in the following lemma. 



Lemma 1. Properties of the fermionic R- operator, 
(i) Evenness. The fermionic R-operator is even, 

p{^^{u,v)) = Q . (103) 

(a) Bilinear relation. The fermionic R-operator satisfies 

^iki^j,^k)^k{u,^k)^j{u,Q = L,{u,^,)C,{u,^k)Oi^.,{^,,^^ ■ (104) 

(Hi) Yang-Baxter equation. The fermionic R-operator satisfies the following 
form of the Yang-Baxter equation, 

:rUu, v)Jl(,{u, w):kUv, w) = Jii,{v, w)Jl(,{u, w):kUu, v) . (lOS) 
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(iv) Regularity. If R{u,v) is regular, say R"^{v,v) = S^5^, then 

^^{v,v) = P,u , (106) 

where Pjk is the graded permutation operator (|J^/ 



(v) Unitarity. If R{u,v) is unitary (of equation (^^), then CRJ^('U, t>) is 
unitary in the sense that 

:Rf^,{u,v)^i^{v,u) = id . (107) 

Proof. (i) The evenness of the fermionic /^-operator is a direct consequence 
of the compatibihty condition (^1]) . 

(ii) Using the commutation relations (|^) and the projection property (|3ll) 
for the matrices Cj^ as well as the compatibility condition ( ^T]) and the 



Yang-Baxter equation (|40|) , the matrix elements on both sides of (|104|) 
can be reduced to 

('_l){p(a)p(/3)+P(/3)p(7)+P{7)p(")+P{/3")(p(7)+P(7"))} 

(iii) The proof is similar to the proof of (ii). Using (|2|), (^ and (|l|), (H) 
both sides of equation ( p.05| ) reduce to 



_-|^Up(a)p(/3)+p(/3)p(7)+p(7)p(a)+p(a"){p(a)+p(a"))+p(/3")(p(7)+p(7"))} 



IV 



(_l)p(7)+p(«)(p(/3)+p(7))5«5^e^.7e^J = (_l)P(/3)e^.^e,^ = P,, 



V 



Using (0), (0) and (|T]) we obtain 



^U^, v)%{v, u) = (-l)Pmp(<^)^PiP'))R-J(u, v)R%,{v, u)ef:e4 
and the assertion follows from (p5|). 



D 
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8 The shift operator 

In this section we shall use the fermionic i?-operator introduced above in 
order to define the shift operator for inhomogeneous graded models. We 
shall explore the properties of the shift operator and shall eventually use 
these properties to prove our main result (pG]). 

Let us start with a slightly more general concept. The inhomogeneous 
monodromy matrix defined in ( p^ is an ordered product of L-matrices. In 
the following we shall indicate the order of the factors by supplying subscripts 
to the monodromy matrix 

7,,„l{u; ei, . . . , ^l) = T(m; ^i, • • • , ^l) = ^l{u, ^l) ■ ■ ■ ^liu, ^^) . (108) 



As can be seen from (|103|) and (|104|) the fermionic _R-operator 'Jijj_^_i{^j, ^j+i) 



interchanges the two neighbouring factors £j+i(u,,^j+i) and Lj{u,^j) in the 
monodromy matrix. Since the symmetric group &^ is generated by the 
transpositions of nearest neighbours, the L-matrices on the right hand side of 
(|108| ) can be arbitrarily reordered by application of an appropriate product 
of fermionic /^-operators. This means that for every r G &^ there exists 
an operator "^l^i^i, ■ ■ ■ ,^l), which is a product of fermionic /^-operators 
and induces the action of the permutation r G &^ on the inhomogeneous 
monodromy matrix, 

The non-graded analogue of this operator was introduced in ||16 |. 



Let us construct "^lii^i, ■ ■ ■ ,C,l) explicitly. We shall use the shorthand 
notation ^R.^ l = ^i l(6, • • • , ^l), 'J'i...l(m) = 'J'i...l(m; ^i, • • • , ^l) and 3lL = 
3?i(0, a) whenever the order of the inhomogeneities 6, • • • , ^ is the same 
as the order of the corresponding lattice sites. For j = 1, . . . , L — 1 define 
aj G &^ by 



a,{k)={j iffc = j + l, (110) 
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The aj G © are transpositions of nearest neighbours. It follows from ( |103|) 
(ITOl) that 



:Kj,+iTi..x(«)=T.,(i)....,(L)(«):J?i+i . (Ill) 

This means that CR^^ ^ = 3?{j+i- Choose r G ©■^ arbitrarily. Then 

Since the transpositions of nearest neighbours cij, j = 1, . . . , L — 1, generate 
the symmetric group 6-^, there is a finite sequence (jp)p=i, such that r = 
Gj^ . . . aj^. Let Tp = (Tj-^ . . .(Tj^, p = 1, . . . , n and tq = id. Then r = r„, and, 
using ( |112| ), we conclude that 



^rj,_iOp),Tj,_i(ij,+l) '^Tj,_i(l) ...rp_i(L)(M) - 'J'tp(I) ...rp(L)(M) ^rj,_iOp),rj,_iOp+l) ' 

(113) 
for p = 1, . . . ,n. By iteration of the latter equation we obtain 

7r(l) r(L){u)yi^ ,., r^U-'-^^f-^ f ^U^^ i +1 ■ (114) 

Thus we have constructed an explicit expression for the operator D^[ j^, 

^l...L = ^T„„iO„),T„_iO„ + l) • • • ^ri(J2),riO-2 + l)^il.il + l • (^-^^) 

Let us now specify to the case, when r is equal to the cyclic permutation 
7 = 0"!... o"L_i. Then jp = p, p = 1, . . . , L — 1, in our above construction, and 
7p = di . . . (7p. Thus 7p-i(jp) = lp~\{v) = 1, and 7p_i(jp + 1) = 7p_i(p+ 1) = 
p + 1. Using ( |115| ) we obtain 

^T...L = ^IL^lL-l • • • ^12 ■ (116) 

The operator "JlJ ^ induces a shift by one site on the inhomogeneous mono- 
dromy matrix. Now ( |112| ) implies that 



'^]{l)...^(L)^^i^)-l{L){.u) = 7^2(i)...^2(i)(M)ai^(^^^^^^(^) . (117) 

It follows by multiplication by "Xl ^ from the right, that 

^l...L = ^7(1)...7(L)^1...L • (118) 

Iterating the above steps we arrive at the following lemma. 
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Lemma 2. The operator 

-^l...L ~ 7"-i(l)...7"-i(L) 7"-2(l)...7"-2(L) • • ' ^l...L ' \^^^) 

where 

•^7P-i(l),..7P-i(L) — -^pp-l ■ ■ ■ -^pl-^pL ■ ■ ■ -^pp+l ' (l^uj 

generates a shift by n sites on the inhomogeneous lattice, i.e. 

^T..L'^l-M = '^n+l...Ll...n{u)^T..L ■ (121) 

Since 7^ = id, we conclude from ( |12lD that 

^t.L ^i...l(^) = 7i...l{u) Jit..L ■ (122) 

If D^jj^. is unitary, we have the following stronger result. 
Lemma 3. Let D^^^ be unitary (cf equation 101 ). Then 



"^L.L = id • (123) 

Proof. Let us first prove the case L = 2. Then D^]^2 = ^12 ^^d CR^j = 

-^7172-^12 — ^21 -'^12 — ^21 -'^12 — id. The last equation holds, since by hy- 
pothesis, 3^(2 is unitary. 



For the case L > 2 we start from the Yang-Baxter equation ( [L05| ), 

^L,L~n'^L,j^L-n,j = '^L-n,j'^L,j'^L,L-n " (124) 

By iterated use of ( |124D we obtain 



■^L,L-n\-^L,L-n-l ' • ' ■^L,l)\-^L-n,L-n-l ' ' ' -*^L-n,l) ~ 

(^L-n,L-n-l- • •^L-n,l)(^L,L-n-l- • •^L,l)^L,L-n ) (125) 



for n = 1, . . . , L — 2. 



Let us introduce the truncated cyclic permutations 7^ = cri . . . (Tp_i, p = 
2, . . . , L, as above. 7^ induces a cyclic shift on the p-tuple (1, . . . ,p) and 
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leaves the [L — p)-tuple {p + 1, . . . ,L) invariant. Using ( |1 25|) , it follows that 

rnf mf <T)7 

^L,L-n ■ ■ ■ -^L,! 7^-1-1(1) ...■y'^-"-'^{L) 

— ■^L,L-n V^L,L-n-l ■ ■ ■ '^L,l)\.-^L-n,L-n-l ' ' ' -^L-n,!) ' 

' \-^L~n,L ■ ■ ■ '^L-n,L-n+l> 

— V^L-n,L-n-l ■ ■ ■ ^L-n,l)v^L,L-n-l ■ ■ ■ -^L,!) ' 

rnf n)f (n)f <T)/ N 

■ •^L,L-n'^L-n,L \'^L-n,L-l ■ ■ ■ -^L-n.L-n+l) 

'• V ' 

= id 

~ \-^L-n,L-n-l ■ ■ ■ •^L-n,l)i-^L-n,L-l ' • • •^L-n,L-n+l) ' 

■ \-^L,L-n-l ■ ■ ■ -^L,!) 
^•^^lzrHi),-,lLZr'iL-l),L-^L,L-n-l---^L,l ■ (126) 



Hence 



L 

— '^ 

7-^-1(1). ..7i-i(L) ■^^7^-2(1). ..7i-2(L) ^1-L 



n?7 _ «T)7 «T)7 «T)7 



mf rnf rp7 rnj nj7 

■^L,L-1 • • • -^L,! •^7L-2(i)_^L-2(2,) -^7^-3(1). ..7i-3(x) • • • •^l...L 

mlL-l rnf rnf fT)7 rT)7 

"fizlm, - ,7^:?(i-l), L '^L,L-2 ■ ■ ■ ^L,l ^7i-3(i)...7i-3(L) ' ' ' ^1. 



(127) 



,L-1 .i-2 



v2 



rn^L-l _ rn^L-2 _ _ fp75 

-^l...L ~ "^l...L — ■ ■ ■ — ^i,,,x, 



Since 72 = (Xi and Di^^ ^ = CR{2) the latter equation reduces the proof of 
Lemma 3 for L > 2 to the case L = 2, which was proved above. D 

Our next lemma can be used to establishes a connection between the inho- 
mogeneous monodromy matrix ( |108| ) and the shift operator ( |119| ) 



Lemma 4. Let X = X^e^ G End(C'""'"") and let R{u,v) be regular, say, 
R?;^iv,v) = 626^. Then 

str(XT„,.xi...n-i(a)) = (-l)^^"^+^(")^^''^X|e„^ai^„_,,,,^„_i^ . (128) 
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Proof. 

Str{X7n...Ll...n-liU)) 

■■■'^n+li3^^^{C,n,C,n+l)y -LJ e„^ 






p(a)+p(«)p(/3)+E-^=i(p(/3,)+p(a,)p(/3,)) 



Y"A^ X/3n-l X/32 Xft r/Sz, X/3n + 2p /3n+l 



p(q)+p{q)p(/3)+E'' ,=i (p{/3j)+p(aj)p(/3j)) 



V" p /3p /3n-lp /?n-2 p Pi p Pl p Pn+2 

'C^n+2/3^^2 (^"' '^"+2) • • • ^Ljs'l {C,n, Cl)^!/?^^ (■Cm Cl) 



p(Q)+p{a)p(/3)+E'' ,=1 (p{/3,)+p{a,)p{/3,)) 
-1)1 j>^»i,"+l 



^+2 

1+2 



. y^ f> f^P Pn~\ p /3n-2 p /Jlp /3i p ^3n-^ 

• '^n+2/3^^2 (^"' ^"+2) • • • '^i/3£ (^«' ^l)'^1/3i ('Cn, Cl) 

• • • ^n~l'^2Zl{^n,^n^l) %i,n+l (129) 

_ (_T\pia)+p{a)p{/3) Y^p P <vf tvf 'T?/ H?/ 

= (-lY('-)+p('^)p(^)X^eJ^3iy,^^^^^^_,^ . (130) 

Here we used the regularity in the first equation. In the second equation we 
reversed the order of factors and introduced a product of Kronecker deltas. 
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In the third equation we used the identity 

"a„_i"Q:„_2 • • • "ai"ai"aL-i ' ' ' "a„+i^"a 

— p /5p Pn-l p /3n-2 p Pi p Pl p Ai + 1 (A'\A\ 

^na^nan-i^nan-2 ' ' ' '^"ai'^^aL ' ' ' ^"Qn+i ' V-^"-"-^/ 

which follows from (|3lD. In the fourth equation we used that 

3^{. = (-ir^'^^+^^-^^^^^e.f £,^(e„ iu) (132) 

and the fact that D^L is even. In the fifth equation we iterated the two 
previous steps of our calculation. Finally in the sixth equation the formula 
(|120|) entered. D 



Setting X^ = (5^ in ( |128| ) and using the cyclic invariance of the super 
trace we obtain the following corollary to Lemma 4. 

Corollary. 

'^yn..n--L = str(Ti...z.(en)) • (133) 

Equation ( |133| ) is the inhomogeneous analogue of equation (|53|) . 



Lemma 5. We have the following expression for the shift operator in terms 
of the elements of the monodromy matrix, 

n 

n:.L = X{M'yi...L{i,)) . (134) 



// i?(M, v) is unitary (of equation 1^ )), then 3^^ ^ is invertible with inverse 



(^i"^)''= n str(0-i...,(e,)) . (135) 

j=n+l 

Proof. The lemma follows from Lemma 2, Lemma 3 and Corrolary 1 to 

Lemma 4. D 



We are now prepared to prove our main result, equation (|86|). 



34 



Proof of equation (8t). Using Lemma 2, Lemma 4, the corollary to Lemma 



4 and Lemma 5 we obtain 
str(X7„...Li...„_i(^„)) 



(136) 



n-l L 

= \{M7i...l{Q) ■ str(XTi...L(en)) ■ ns*^(^i-^(0)) 

j=l j=n 

= (-l)P(")+^(")P(^)X|e„^str(Ti...z.(a)) ■ 
It follows that 



n-l L 

nstr(Ti...L(0)) ■str(XTi...L(en))- J] str(Ti...i(0)) • (137) 

j=l j=n+l 

Finally, by specifying X^', = (— 1)p(°')+p(°')p(/3') ^^'5^,, we arrive at equa- 
tion (m. n 



Summary 

In this article we obtained an explicit solution of the quantum inverse problem 
for fundamental graded models. Our main result is the general formula (^61). 
This formula expresses the local projection operators e„(^, which represent 
local spins and local fields, in terms of the elements of the monodromy matrix. 
The formula and its proof essentially simplify for translationally invariant 
models (all inhomogeneities coincide, C,j = 0- ^^ ^^e translationally invariant 
case the proof is based on the representation of the shift operator as a product 
of permutation matrices. 

We presented explicit formulae for the solution of the quantum inverse 
problem for the XYZ quantum spin chain (p!7D-(p!9|) and for the supersymmet- 
ric t-J model of strongly correlated electrons (|9ll)-(^8D. For the XYZ chain 
the local projection operators coincide with local Pauli matrices (quantum 
spin operators). 
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We are planing to use our results in forthcoming publications in order 
to obtain multiple integral representations for correlation functions of funda- 
mental graded models. 
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